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ABSTRACT 

We present the observed offsets of short-duration gamma-ray bursts (SGRBs) from 
their putative host galaxies and compare them to the expected distributions of merging 
compact object binaries, given the observed properties of the hosts. We find that for all 
but one burst in our sample the offsets are consistent with this model. For the case of 
bursts with massive elliptical host galaxies, the circular velocities of the hosts' haloes 
exceed the natal velocities of almost all our compact object binaries. Hence the extents 
of the predicted offset distributions for elliptical galaxies are determined largely by 
their spatial extents. In contrast, for spiral hosts the galactic rotation velocities are 
smaller than typical binary natal velocities and the predicted burst offset distributions 
are more extended than the galaxies. 

One SGRB, 060502B, apparently has a large radial offset that is inconsistent with 
an origin in a merging galactic compact binary. Although it is plausible that the host 
of GRB 060502B is mis-identified, our results show that the large offset is compatible 
with a scenario where at least a few per cent of SGRBs are created by the merger of 
compact binaries that form dynamically in globular clusters. 
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1 INTRODUCTION 

The merger of a compact binary consisting of neutron stars 
(NSs) or black holes (BHs) is a prime model for the origin of 
the short-duration gamma-ray bursts (SGRBs). Compact bi- 
naries may be formed by numerous channels includi ng those 
which are formed in primordial binaries (Bclcz vnski et al.l 
|2002| ). and those which are formed dynamically in dense 
cluster cores (e.g. lDavieslll995l : ICrindlay et ai] 120061 '). Once 
formed the binary loses energy and angular momentum from 
its orbit in the form of gravitational radiation, and the 
compact objects spiral together and eventually merge. The 
merger time is highly sensitive to both separation and ec- 
centricity and thus has an extremely broad distribution. 

Evidence for the origin of SGRBs in the final co- 
alescence of such systems potentially comes both from 
the host galaxy types (e.g. IZheng fc Ramirez-Ruij |2007| ) 
and the measured offsets of SGRBs from their hosts (e 



iBloom et alll 19991 : iBelczvnski et al.ll2006l : iTroia et al.ll2008 

The kicks imparted to double compact object (DCO) sys- 
tems on formation, both in terms of dynamical kicks to the 
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binary and natal kicks to neutron stars, endow DCO binaries 
with spatial velocities of up to several hundred kms~^. This 
means that they may merge far from their birth sites, even 
outside their parent galaxies. Equally, if the inspiral takes 
> 10* years, then the host galaxy may no longer exhibit sig- 
nificant ongoing star formation at the time of merger. 

Observations of SGRBs obtained by Swift so far of- 
fer broad support for t his picture. Several of the host 
galaxies observed are old (ICehrels et a l. 2005; Berger et al.l 
l2005l ; iBarthelmv et al.ll2005l ; iBloom et al.i i2006). while only 
a small fraction are starbursts. However, the sample has 
been building up slowly, and we are only now approach- 
ing a stage where it is possible to compare the observed 
offsets with predictions for DCO binaries. We build on pre- 
vious work tha t model led the offset distributions of DCO 
mergers (e.g. iBloom et al . 1999; Fryer ot al. 199jJ) by using 
the observed properties of the host galaxies to predict offset 
distributions on a host-by-host basis. 

In this paper we take the observed sample of SGRBs 
that have identified hosts and model the production and 
galactic trajectories of DCO binaries inside those hosts. In 
Section 2 we review the sample of bursts; we describe our 
models in Section 3, discussing the results and possible se- 
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Table 1. Properties of SGRBs with known X-ray, optical or radio counterparts. Data is taken from the literature sources indicated 
other than for a small number of the offsets which we have re-measured in described in Section 2. Columns list in order: the burst 
identifier, the satellite which detected the burst, the duration over which 90% of the total fluence was seen (tgo), whether an optical 
counterpart was detected, the host apparent i?-band magnitude and error, the host absolute B-band magnitude, the redshift, the host 
type (spiral/elliptical), the host effective radius Rc, the offset fiproj and the error on the offset Rcrr- The duration, igo, is measured in the 
15-350 keV range for the Swift bursts and in the 30-400 keV range for the HETE-2 burst. The absolute magnitudes have been calculated 
assuming a source with a flat spectrum in . All offset errors are o ne-sigm a; 90% confidence limits in the literature have been converted 
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lection effects in Section 4. The alternative possibility of 
DCO binaries forming dynamically in the cores of globular 
clusters is discussed in Section 5. 



2 THE BURSTS AND THEIR HOST 
GALAXIES 

From the set of observed short GRBs we have selected the 
bursts with positions localised through observations of their 
X-ray and, where possible, optical afterglows. Furthermore 
we have discarded bursts without identified host galaxies or 
for which there are not robust redshift measurements, since 
knowledge of the redshift is necessary to construct model 
galactic potentials. For several bursts we have re-analysed 
available data to perform relative astrometry at times when 
the afterglow was bright, mapping the locations of the burst 
onto the hosts utilising the geomap task within IRAF and a 
number of point sources within the field. For other bursts we 
have used offsets available in the literature. Our final sample 
of 16 bursts is listed in Table [1] 

In order to obtain dynamical models for the bursts' 
host galaxies we split the sample into elliptical and spiral 
hosts. In both cases we adopt the the logarithmic profile of 
omas et al.l (|2009| ) to model the host galaxy's dark halo. 



vi, 3rg + 
47rG {fl -f r-2) 



(1) 



where rh is the core radius of the halo and its circular 
velocity at infinity. This profile is constructed such that the 
circular velocity approaches Uh asymptotically at infinity. 



It should be noted that as elliptical galaxies do not have 
discs in their case the circular velocity is not directly mea- 
surable, but it is useful to parameterise the potential in 
the same way. Initially we place the binaries in an expo- 
nential disc, chosen so that half the stars are within the 
effectiv e radius. For el l iptica l hosts we use the scaling rela- 
tions of lThomas et al.1 (|2009l ) to obtain the properties of the 
galaxy's halo from its bl ue-band magni t ude. The effective 
radius Re is taken f rom iGerhard et al.l (|200ll ) as cited by 
iThomas et alj ll2009h. The proper t ies of spiral galaxies are 
taken from Kormendv fc FreemanI (|2004l ). The complete set 
of models that we used are given in Table[21 We take f?e to be 
the half-light radius of the galaxy. A comparison between the 
observed and predicted Rc shows that there is considerable 
scatter, as would be expected from observations o f galaxy 
properties (see e.g. figure 3 of lThomas et al.ll2009l ). Where 
a galaxy has an observed Re we compute models using both 
the value predicted by the formula and the measured value, 
and use the observational value in preference (for example, 
in the cumulative offset distributions of Figures [SI and [TU)) . 



3 POPULATION SYNTHESIS AND 
GALACTIC TRAJECTORIES 

To model the observed DCO populations we conduct a popu- 
lation synthesis of NS-NS and BH-NS binaries , utilising the 
rapid binary population synthesis code B SE jHurlev et al 



I2002iy We have modified B S E fol lowing iBelczvnski et al 



1 20021 ) and IBelczvnski et al.l (|2008l ) to include a more re- 
alistic prescription for compact object masses, hypercritical 
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We specify the natal kick distributions for the two ob- 
jects that form in each DCO separately, considering three 
combinations of kick distributions. In the first scenario both 
compact objects receive an ACC02 (strong) kick. Alter- 
natively, the first compact object to form may receive an 
ACC02 (strong) kick and the second a DPP05 (weak) kick. 
Finally both objects may receive a DPP05 (weak) kick. As 
the evidence for weak kicks relies on the properties of the 
second compact object we do not consider the case where 
the first kick is weak and the second kick is strong. 

A powerful test of the neutron star kick distributions 
is the comparison of the synthesised NS-NS binary popu- 
lations that they yield with the properties of observed NS- 
NS binaries. In particular, the kicks very strongly affect the 
semi-major axis a and eccentricity e of a given DCO binary. 
The ability of a kick distribution to produce the observed 
NS-NS binaries in the a, e plane is a stringent test of its cor- 
rectness. We present the a, e distribution of all the NS-NS 
binaries formed in our simulations in Figure [T] 

There is little evidence as to whether a kick is im - 
parted to black holes at formation (|Belczvnski et al]|2008l '). 
though for the Galactic black-hole binary Cygnus X-1 
there is evidence that there was no substantial natal kick 
(|Mirabel fc Rodrigued l2003l) . A significant kick would re- 
quire the asymmetric ejection of a large fraction of the pro- 
genitor's envelope, rather than its accretion onto the black 
hole. The larger mass of a typical black hole compared to 
a neutron star increases the quantity of mass that must be 
ejected to produce a significant kick. Hence we choose not 
to apply a natal kick to black holes upon formation. Where 
a black hole forms via the accretion-induced collapse of a 
neutron star we of course apply a kick to that neutron star 
at its birth. We have tested the effect of omitting black hole 
kicks by synthesising a population of BH-NS binaries where 
a kick following the ACC02 distribution is applied, reduced 
according to the ratio of the typical neutron star mass to 
the nascent black hole mass. This is equivalent to a natal 
kick with the same distribution of impulses as a neutron 
star kick. We found that the space velocity distribution of 
the binaries that we produced was largely unaffected. This 
is unsurprising as such binaries are typically rather massive 
and hence strongly bound at the time of the kick and, as our 
results show, the properties of our sample of DCO binaries 
is relatively insensitive to the natal kick chosen for the first 
supernova. 

In addition to the natal kicks each supernova causes an 
impulse to the binary by the rapid loss of a large quantity 
of mass. As this mass escapes anisotropically it takes away 
a significant net momentum. This effect is also taken into 
account in our simulations and accounts for the fact that, 
even when two small kicks are utilised, a DCO binary may 
have a large spatial velocity after formation. 

We plot the mass lost in the final supernova as a func- 
tion of the orbital separation immediately prior to the super- 
nova for all potential DCO progenitors in Figures [2] (NS-NS 
binaries) and [3] (BH-NS binaries). In each case the upper 
panel shows binaries that are bound after the supernova 
and natal kick, whilst the lower panel includes systems that 
unbind at the final supernova and hence form two isolated 
compact objects. In both cases the natal kick distribution is 
that which we prefer for the subsequent analysis, the ACC02 
distribution (see Section [4.11) . It is apparent in both cases 



Table 2. Parameters of the galactic models that we used. 
Columns list the GRB, whether the host is spiral or elliptical, 
the asymptotic rotational velocity of the halo model vjj, the halo 
core radius r^, the predicted effective radius Rc and, where one is 
available, the observed effective radius R'^^ . The final column is 
the mass enclosed within a spherical radius of 10 -Rc, which acts 
an an indicative halo mass. 

accretion during common envelope evolution, and delayed 
dynamical instability in mass transfer from helium stars in 
binaries with a large mass ratio. More details of our popula- 
tion synthesis code, and the constraints which can be placed 
on DCO formation will be presented in a later paper (Church 
et al. in prep.). For reas ons of simplicity we d istribute stellar 
masses according to the lKroupa et all (|l993t l IMF, selecting 
both stars independently from the IMF and only consider- 
ing stars with masses greater than 3 Mq . Binary semi- major 
axes are chosen from a distribution fiat in log a between 1 
and lO^i?0. Preliminary simulations showed that binaries 
wider than this do not contribute to the DCO population. 
The initial eccentricity of the binaries has little effect on 
their evolution and is set equal to 0.1 for all binaries. 

3.1 Neutron star kicks 

For the purpose of modelling the galactic offsets of SGRBs 
the most significant uncertainty is the natal kick imparted 
to neutron stars at their formation. In this paper we utilise 
two distributions of neutron star natal kicks: the bimodal 
kick distribution of lArzoumanian et all l|2002l ). henceforth 
ACC02, and the Maxwellian distribution recommended by 
iDewi et all (|2005l '). henceforth DPP05. 

The ACC02 distribution leads to strong kicks, with typ- 
ical magnitudes of several hundred kms~^, and is derived 
from observations of isolated pulsars. The DPP05 distri- 
bution yields much smaller kicks - the dispersion is only 
20kms~^ - and is derived from constraints placed on the 
space velocity of the second neutron star formed in a binary 
by the relationship between its spin and orbital eccentric- 
it y. We also computed mode ls that utilised the distribution 
ofH ansen k, PhinnevI (| 19971 ). but this produced results that 
were, for our purposes, very similar to those obtained with 
the ACC02 distribution. 
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Figure 1. Dots show the eccentricity e as a function of semi- 
major axis a for our synthesised populations of NS-NS binaries, 
plotted in the location in which they form. To clarify the plots 
only 2000 binaries are plotted in each case. Crosses show the 
properties of the observed NS— NS binaries. Lines of iso-merger- 
time of 10*^, 10* and 10^° yr are also plotted. The three panels 
show different combinations of kick distributions. Top: ACC02 
kick for both NSs. Middle: ACC02 kick for the first NS and DPP05 
kick for the second NS. Bottom: DPP05 kick for both NSs. 



that the pre-supernova binaries fall into two groups, one with 
large separations and one with small separations. The large 
separation group tend to also have larger mass-loss, as they 
have not gone through a recent episode of common- envelope 
evolution and hence the progenitor of the second compact 
object retains a larger fraction of its envelope. These wider 
binaries are usually disrupted by the final supernova, un- 
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Figure 2. Dots show the mass lost in the final supernova, AM, 
as a function of pre-supernova semi-major axis a for the binaries 
that potentially form a NS— NS binary. The bottom panel shows 
all binaries that are bound prior to the final supernova, the upper 
panel only the subset that remain bound subsequent to the final 
supernova. The ACC02 kick distribution was used for both kicks. 

less a chance alignment of natal kick causes them to remain 
bound. In the case where they do produce a bound binary it 
is usually wide enough that it does not merge under the ef- 
fects of gravitational wave radiation within the Hubble tim^B 
and hence is not significant for our calculation. The major- 
ity of systems that go on to merge and potentially form a 
gamma-ray burst come from the population of closer bina- 
ries consisting of a compact object and a helium star. 

The merger time distributions for the three different 
models are shown in Figure |4l Under all three assumptions 
about the strength of kicks, the merger times are predom- 
inantly long; that is, greater than 1 Myr. The distributions 
that we obtain are similar to distributio ns obtained by p re- 
vious work in the field (e.g. figure 23 of lFrver et al.|[l999l ). 

3.2 Galactic trajectories 

In order to predict the offsets of compact binary systems 
from their host galaxies at the time of final merger we in- 
tegrate their motion in the gravitational potentials of their 

^ The Hubble time is taken to be 14Gyr throughout. 
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Figure 3. As in Figure[2]but for BH-NS progenitors. Dots show 
the mass lost in the final supernova AM as a function of pre- 
supernova semi-major axis a for the binaries that potentially form 
a NS-NS binary. The bottom panel shows all binaries that are 
bound prior to the final supernova, the upper panel only the sub- 
set that remain bound subsequent to the final supernova. The 
ACC02 kick distribution was used for both kicks. 




modelled host galaxies. For each host we produce a halo 
model as described in Section (2] the parameters of the mod- 
els are listed in Table [2] We place the stars initially in an 
exponential disc and assign the initial location of each DCO 
binary with a probability proportional to the disc's lumi- 
nosity at that point, i.e. our DCO binaries originate tracing 
the light of their host galaxies as is seen for most supernovae 
l|Fruchter et~aIll2006l : lKellv et al.li200 8'). Each DCO binary's 
initial velocity is set equal to the sum of its circular velocity 
in the galactic potential and a kick velocity that originates in 
the effects of the supernovae on the binary's centre-of-mass 
velocity; this is calculated during the population synthesis. 
The orientation of the orbital plane of the binary with re- 
spect to the disc of the host galaxy is taken to be isotropi- 
cally distributed. The binary's trajectory is then calculated 
and its position at the time of merger recorded. We view the 
galaxy from a direction chosen uniformly across a sphere and 
record the binary's projected offset from the galaxy centre 
at this time. We compute 10^ trajectories for each galactic 
model in order to obtain distributions of predicted burst off- 



I .^■'—r-^ -T , T- - , I . , . lLJ 

3 4 5 6 7 8 9 10 

logio 7-merge/yr 

Figure 4. The distribution of merger times of the NS— NS and 
BH-NS binaries produced by our population synthesis calcula- 
tions. Solid lines show NS-NS binaries, dotted lines BH-NS bi- 
naries. In the top panel the ACC02 kick distribution was used 
for both kicks, in the middle panel the ACC02 distribution was 
used for the first kick and the DPP05 distribution for the second 
kick, and in the bottom panel the DPP05 distribution was used 
for both kicks. Only binaries that merge within the Hubble time 
are included. 



sets. These distributions are shown in Figures (5] [6] and [71 in 
addition to the offsets of observed bursts. 

For some of the hosts the inclination of the host disc 
with respect to our galaxy is known. Hence we tested 
whether the viewing the host from a fixed angle as opposed 
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to a random angle affected our results. We found that the 
distribution of offsets did not change significantly. 
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4 DISCUSSION 

The results of the simulations are presented in Figure [Tl 
which shows the distribution of binaries in the a vs. e plane, 
and Figures [5] to [T] which contain the binary natal velocity 
distributions and predicted offset distributions. 



4.1 Synthesised binary population 

The distributions of compact object binaries in the a, e 
plane that our population synthesis calculations generate are 
shown in Figure [T] A detailed discussion of whether the ob- 
served population is consistent with our prediction requires 
consideration of the biases in the detection of pulsars and 
is beyond the scope of this project, but it can be simply 
stated that it is necessary for all of the observed binaries 
to be able to form under the adopted description of binary 
evolution processes. The a, e plots show that, for a kick dis- 
tribution where the second kick is weak - that is, drawn 
from a DPP05-like population - it is not possible to make 
a NS-NS binary such as PSR B1913-f 16 (the Hulse-Taylor 
pulsar). This means that the binary in which it formed must 
have experienced a relatively strong second kick. It does not 
necessarily indicate that all NS natal kicks must be strong, 
as some intricacy of the evolution may cause some NSs to 
receive a strong natal kick and some a weak one, but for the 
purposes of this study it leads us to prefer a kick distribu- 
tion where all the NSs receive a strong kick on formation. 
The a-e distribution that we obtain using a stron g kick is 
similar to that of, for example. iBloom et al.l (| 19991 ). though 
with an additional population of long-period binaries. This 
difference is not significant for this project as the long-period 
population all have merger times greater than the Hubble 
time. 

Examination of the binary birth velocity distributions 
presented in Figures[5]to[71 which are the result of the effects 
of both supernova mass loss and supernova natal kicks on the 
space velocity of the compact binaries, shows that there are 
significant differences between NS-NS and BH-NS binaries. 
In the cases where the second kick is weak (Figures[S]and[71) 
a significant population of the BH-NS binaries acquire ve- 
locities that are large relative to the NS natal kick velocities. 
These large velocities originate in the loss of a substantial 
quantity of mass during the formation of the NS when the bi- 
nary is close and hence the stars' velocities large. The lower 
compact companion mass in comparable systems that form 
NS-NS DCOs limits the amount of mass that can be lost 
without the binary breaking up and hence reduces the ef- 
fect of this process. Therefore, for a small second kick the 
systems with the largest space velocities are BH-NS bina- 
ries. This is not the case, however, when the second kick 
is large (and hence more dominant over the effect of mass 
loss); there a much larger proportion of the highest-velocity 
systems are NS-NS binaries. 
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Figure 5. Top: the distributions of birth velocities for NS-NS 
(solid lines) and BH-NS (dashed lines) binaries. Both first and 
second natal kicks are taken from the ACC02 distribution. Mid- 
dle: offset distributions of NS-NS binaries around different burst 
host galaxies compared to the observed burst offsets. The dark- 
ness of the colour is linearly related to dP(i?)/dlogi?. Error bars 
show one-sigma uncertainties in the observed offsets. Where both 
observed and predicted effective radii are available the additional 
right-hand stripe shows calculations made with the observed iJo . 
The two rightmost stripes represent the distribution predicted for 
the globular cluster systems of M87 and the Milky Way. Bottom: 
as middle but for BH-NS binaries. 
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Figure 6. As in Figure [5] but the first l<icl< is taken from the 
ACC02 distribution and the second from the DPP05 distribution. 



Figure 7. As in Figure \5\ but with both kicks drawn from the 
DPP05 distribution. 



4.2 Predicted galactic offsets 

The hosts in Figures [S] to [7] are sorted in order of decreas- 
ing blue-band magnitude. A generally large scatter can be 
seen in their observed offsets, though two of the three most 
luminous - and more massive - elliptical galaxies host short 
bursts at typically large offsets. These are also typically 
bursts with shorter durations, which have been suggested 



to lie at larger radii from their hosts a et al.l |2008| ). 

where the short durations, and very faint X-ray afterglows 
would correspond to low ISM density about the progenitor. 
The large scatter observed is consistent with the models, as 
shown by the broad probability distributions. 

There is a qualitative difference between the behaviour 
of binaries in the elliptical galaxies in our sample and those 
in spiral galaxies. The elliptical galaxies are much more 
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Figure 8. The probability of obtaining an offset larger than some 
projected radius r as a function of r, in our models of the host 
galaxy of GRB 060502B. Solid lines represent the distributions 
of NS-NS binaries, dashed lines those of BH— NS binaries. Distri- 
butions have been calculated using the observed _Re- The three 
panels show the three different pairs of kick distributions that we 
considered. The error bars, placed at arbitrary heights, show the 
one-sigma error on the observed offset. The dotted line marks the 
three-sigma probability. 



massive and have halo circular velocities between 400 and 
700kms~^. These circular velocities exceed the vast major- 
ity of the binary natal velocities regardless of the kick dis- 
tribution used. Hence the binaries' galactic orbits are only 
mildly perturbed by the kicks and the offset distributions 
are determined largely by the effective radii of the hosts. 
The spiral galaxies, on the other hand, have much lower cir- 
cular velocities of between 100 and 200kms~^. Hence they 
have extended offset distributions given their small sizes, 
particularly for the BH-NS binaries and the stronger kick 
distribution. 

Our models predict offset distributions that match all 
of the bursts observed around spiral galaxies, provided that 
either the natal kick is relatively strong, or the progenitors 
of the bursts with the largest offsets are BH-NS binaries 
rather than NS-NS binaries. The cumulative offset distribu- 
tions for the spiral galaxies hosting the two bursts with the 
largest offsets are shown in Figures [g] (GRB 061217) andfUJl 
(GRB 061201). However, the burst with the largest offset, 
GRB 060502B which has an elliptical host, is not well fitted 
by the models. The cumulative burst probability, calculated 
using the measured -Rc , is shown in Figure [S] Fewer than 
0.1% of possible hosts lie within its 1-sigma offset position, 
and fewer than 2% within its 3-sigma offset position. We 
cannot exclude the model at a 5-sigma level as the error box 
then overlaps a large fraction of the host galaxy. 

Given the uncertainty in the distribution of natal kicks 
it is prudent to investigate the effects of a natal kick stronger 
than the ACC02 distribution. Following a suggestion from 
the referee, we re-calculated the cumulative offset distribu- 



Figure 9. As Figure |8] but for GRB 061217. The calculations 
use a predicted Re as no observation was available. 
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Figure 10. As Figured but for GRB 061201. 



tion using an ACC kick distribution scaled up by a factor 
of two - i.e. where the peak probability was at twice its 
normal value. The results are presented in Figure [11] The 
distribution of offsets from NS-NS mergers remains effec- 
tively unchanged as the stronger kicks merely serve to break 
up a larger fraction of the binaries. A small tail of BH-NS 
mergers at larger offsets is generated, but it is not significant 
enough to explain the burst. 

Several other bursts also lie towards the edge of the pre- 
dicted distribution. Only 8% of predicted burst locations lie 
outside the measured offset of GRB 061201, and only 16% 
outside the measured offset of GRB 070809 (although in 
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Figure 11. As Figure |8] but showing the effect of doubUng the 
strength of the ACC02 kick distribution. The plot shows the 
probability of obtaining an offset larger than some projected ra- 
dius r as a function of r, in our models of the host galaxy of 
GRB 060502B. Solid lines represent the distributions of NS-NS 
binaries, dashed lines those of BH-NS binaries. Distributions have 
been calculated using the observed Re . The three panels show the 
three different pairs of kick distributions that we considered. The 
error bars, placed at arbitrary heights, show the one-sigma error 
on the observed offset. The dotted line marks the three-sigma 
probability. 

that case the positional error bars are larger). Furthermore, 
there is some u ncertainty ov e r the properties of the host of 
GRB 050509B. iBloom et aP (|2006l ) give its effe ctive radius 
as 3.4 7 kpc, rather than the 20.98 kpc given by iFong et al.l 
(|2009l ). This, in combinat ion with the velocity dispersion, 
which iBloom et all (|2006l ') measure to be 250 km s , sug- 
gests that it may be more like the host of GRB 060502B 
than the model we have used. If so this would provide an- 
other large offset that would be hard to match with the 
DCO model. These points strengthen the case that there is 
an excess of large offsets compared to our predictions. 




Projected radius r 



Figure 12. The cumulative probability distribution of projected 
burst offsets plotted for different times after a starburst. Solid 
lines represent bursts occurring within 100 Myr, dashed lines 
bursts occurring between 100 Myr and 1 Gyr after the starburst, 
dotted lines bursts occurring between 1 Gyr and 10 Gyr after the 
starburst and dash-dotted lines burst occurring more than 10 Gyr 
after the starburst. The upper panel shows the distributions for 
NS— NS progenitors, the lower panel those for BH-NS progenitors. 
The model used is that for the host of GRB 060502B. 



4.3 Selection effects and correlations 

There are clearly possible selection effects in operation. 
For example, bursts expelled large distances from low mass 
galaxies may be essentially impossible to associate with 
their hosts with any confidence. This is the likely explana- 
tion for the appare ntly host-less GRBs 061201 and 080503 
jPerlev et al.ll200^ ). A further selection effect worthy of in- 
spection is the possibility that the majority of DCOs that 
form in elliptical galaxies have merged before we observe 
their hosts. If the DCOs that merge at small offsets come 
predominantly from the population with short merger times 
then we would not expect to see them in nearby, old giant 
ellipticals. To investigate this we have plotted the distribu- 
tion of projected burst offsets for our model of the host of 
GRB 060502B as a function of time since a single star forma- 
tion burst (Figure [T^ . The distributions for different times 
are essentially identical, as the galaxy rotation velocities are 
much larger than the supernova kick velocities. Furthermore, 
the fraction of SGRBs observationally associated with old, 
rather than young populations is at least ~ 25%, based on 
the fraction of identified hosts that are elliptical galaxies. 
This strongly suggests that the rapidly merging population 
cannot dominate at the current epoch. 

Other correlations of burst properties with those of 
their hosts have been reported elsewhere. In particular, 
iTroia et al.l l|2008f ) suggest that bursts with extended emis- 
sion typically form at smaller projected distances from their 



hosts than very short bursts. They suggest that this may 
be due to the bursts with extended emission occurring via 
NS-BH mergers, which may travel less far. Our population 
synthesis indicates that this is not the case, and that BH-NS 
and NS-NS populations have broadly similar radial distri- 
butions, with the BH-NS population being more extended 
as explained in Section 14.11 For the preferred strong kick 
we find that the distributions are almost identical (see Fig- 
ure [5]) . 

A second reported correlation was that of 'Rho adsl 

(|2008h . who showed that there is an anticorrelation between 
the absolute magnitude of the host, and the isotropic en- 
ergy release of the burst. Our results confirm this, however, 
in the light of previous discussion it is possible to recast 
this correlation in terms of the observed offsets from host 
galaxies, notably, that the offset and energy appear to be 
anti-correlated as well. It is therefore interesting to inves- 
tigate if the source of the correlation is the due physically 
to the difference in host absolute magnitude, the offset, or 
another physical mechanism which can explain both. 

There is a correlation between the direction of the natal 
velocity of a DCO binary and the orientation of its orbital 
plane. In particular, if the natal kick is weak and hence the 
binary natal velocity is dominated by the loss of mass in the 
second supernova explosion, then the natal velocity should 
be perpendicular to the orbital axis. If the gamma-ray emis- 
sion is beamed along the orbital axis then this correlation 
has the potential to increase the magnitude of the observed 
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offsets. This effect was considered at an early stage in these 
calculations and found to be small once realistic binary evo- 
lution and galactic trajectories were taken into account. 

4.4 Host galaxy selection 

An ongoing problem with studies of SGRBs is the associ- 
ation of a given burst with its host. In a few cases opti- 
cal afterglow positions provide high confidence associations 
showing that a burst lies on the stellar field of a given 
host. However, in many cases this was not possible, and 
different approaches h ave been tak e n. Fo r example, in the 
case of GRB 060502B, iBloom et al] (|2007l ) identify the host 
as the giant elliptical con sidered here with a large offset, 
whilst iBerger et al.l (|2007l ) adopt a strategy of identifying 
the brightest extended object withm the XRT error box, 
which leads to a different, fainter host at a smaller off- 
set. This problem is most severe for GRBs 050509B and 
060502B. The probability of random association with a 
background galaxy some distance from the burst position 
is frequently rather similar to t he probability of a chance 
alignment within the error circle l|Levan et al.H2007l ). In such 
circumstances it becomes very difficult to identify the host 
with high confidence, and it is likely that strong constraints 
will only come with the build up of a larger sample, in which 
it is possible to say if, for example, large elliptical galaxies 
are over-represented. 

Of course, in part the difficulty in locating burst after- 
glows may be due to the progenitors of the bursts them- 
selves. If NS-NS binaries are expelled from their hosts with 
large velocities then their optical afterglow brightness will be 
significantly suppressed due to the low ISM density, since the 
afterglows are caused by external shocks. In very low den- 
sity media the cooling break may even move abo ve the X-ray 
band, and result in very faint afterglows. Indeed. pTroi a et al.l 
(|2008l ) have suggested that bursts with shorter durations, 
and fainter X-ray afterglows lie systematically at larger dis- 
tances from their hosts than the longer duration brighter 
bursts. 

A recent and closely-related paper jBergeJlioiol ) has 
discussed the possible origin of short GRBs with optical af- 
terglows with significant offsets. They concentrate on bursts 
with optical afterglows as having smaller positional errors, 
and show that for all but one burst where the identity of 
the host is in doubt the lowest probabilities of chance co- 
incidence is associated with bright galaxies at offsets of a 
few tens of kpc. Their results are in agreement with ours, 
in that they state that the observations fit a model where 
these offsets arise naturally in merging NS-NS binaries that 
have been kicked out of their host galaxies. They do not 
find any additional contribution at larger offsets necessary; 
on the other hand they do not carry out host-by-host mod- 
elling as we do and they neglect the bursts with only X-ray 
detections that have the largest offsets. 



5 GLOBULAR CLUSTERS 

An alternative explanation of the very large offsets of some 
SGRBs from their host gala^xies is that they occur in com- 
pact binaries residing within the globular cluster systems of 
their host galaxies. Within the dense environments at the 
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Figure 13. The long-dashed thick line shows the offset distribu- 
tion of SGRBs, calculated under the assumption that they occur 
via the merger of dynamically-generated NS— NS binaries in the 
cores of globular clusters and hence follow their projected spa- 
tial distribution. The globular cluster population chosen is that 
of M87. The error bar, placed at arbitrary height, represents the 
offset of GRB 060502B from its putative massive elliptical host 
galaxy. The thin solid and short-dashed lines are the results from 
the standard model of a field population of DCO binaries as pre- 
sented in FigureO The thin solid line shows NS-NS binaries and 
the thin dashed line BH— NS binaries. 

cores of globular cluster s NS-NS binaries can form through 
dynamical interactions jPaviej [l995). Such a scenario for 
t he fo rmation of SGRBs was proposed by iGrindlav et al.l 
(:200(|). who show that, based on some simple assumptions 
about the formation channels of NS-NS binaries within GCs 
and plausible number densities of neutron stars and low- 
mass X-ray binaries, a significant fraction of SGRBs could 
form via this channel. 

5.1 Radial distribution 

To assess the viability of this alternative, we have plotted 
the galactic distribution of gamma-ray bursts assuming that 
they follow the globular cluster population of two galaxies; 
M87 and the Milky Way. M87 is a giant elliptical galaxy 
in the Virgo clus ter, with a population of roughly 14000 
globular clusters l|Harrisll2009l ). The distribution of its clus- 
ters has a very large spatial extent, extending to at least 
100 kpc. The Milky Way has a more modest distribution 
of about 160 clu sters with the most distant also being at 
roughly 100 kpc l|Diorgovski fc Mevlanlll994 ). We plot the 
projected radial distributions of the two globular cluster sys- 
tems in Figu res [S] to [7] For M87 we utilise the fits given by 
iHarrid (|2009l '). whereas for the Milky Way we take the fit 
to the 3-dimensional number d ensity distribution provided 
bv lDiorgovski fc MevlanI (| 19941 ) and project it. For simplic- 
ity we assume that the dynamical formation of the NS-NS 
binaries which may lead to SGRBs proceeds in a similar 
manner across all the GCs. Hence we obtain projected off- 
set distributions, which are plotted on the right-hand side 
of Figures [5] to [T] along with the cumulative probability dis- 
tribution in Figure 1131 Both distributions are much more 
extended than the kick-derived offset distributions, and are 
consistent with the most-offset bursts. 

5.2 Formation rates 

To assess the viability of globular cluster systems as for- 
mation sites for SGRB progenitors it is instructive to esti- 
mate the formation rates within the clusters and galaxies. 
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To simulate the populations here we evolved 4 x 10^ bina- 
ries, which when using strong (ACC02) kicks produced 1659 
NS-NS binaries and 4253 BH-NS binaries, counting only 
those binaries that merge within 1.4 x 10^° years. Once we 
have corrected for the fact that we only sample the heavy 
end of the IMF this leads to specific merger rates of 15 /b 
and 38 /b (10® Mq)~^ ^uubbie- Here /b is the fraction of stel- 
lar mass in binary systems. For globular clusters we extrap- 
olate from the observation of a single dynamically-formed 
NS-NS binary in a Galactic globular cluster, M15-C, with 
a merger time of roughly 300 Myr, amongst the roughly 
150 Galactic globu lar clusters. Utilising the observations of 
iRhode et al.l l|2005h which suggest that the number of glob- 
ular clusters per 10® Mq of stellar mass, T, lies between 1 
(for field spiral galaxies) and 4 (for cluster elliptical galax- 
ies) gives a NS-NS merger rate from globular clusters of 
O.3ir(lO®M0)"^fHibbio- This is obviously much smaller 
than the field merger rate even for elliptical galaxies but 
there is a severe underestimate for several reasons. Pulsar 
surveys of galactic globular clusters are incomplete, and the 
inclusion of only a single system in our calculations sam- 
ples only a small range of possible merger times. Based on 
sim ple scattering simulati ons of binaries in globular clus- 
ters iGrindlav et al.l (|2006t ) show that this simple extrapo- 
lation is an underestimate by a factor of between 10 and 
100, which brings the rate estimate to a similar magnitude 
to that from field binaries. Furthermore, the massive ellip- 
tical galaxy M87 has a va lue of T ~ 8 for the bl ue globu- 
lar cluster sequence alone (|Brodie fc StradeillioO^ ). Clearly 
more detailed simulations of globular clusters are required 
to form a more accurate estimate of their production rates of 
NS-NS binaries, but this crude analysis shows that the rates 
could be comparable, and that given a sample of nearly 20 
SGRBs it would not be surprising for the progenitor binary 
of one or more of them to have formed in a globular clus- 
ter. Furthermore, the burst that we suggest is most likely to 
have formed in a globular cluster, GRB 060502B, has a gi- 
ant elliptical host which would be expected to have a larger 
specific globular cluster fraction. 



6 SUMMARY 

We have synthesised populations of NS-NS and BH-NS bi- 
naries and used them to predict the galactic offset distribu- 
tions of short gamma-ray bursts, given the observed prop- 
erties of their host galaxies. The offsets of all but one short 
gamma-ray burst are found to be consistent with this pic- 
ture. However, the burst GRB 060502B apparently has an 
inconsistently large offset. It is plausible that this inconsis- 
tency results from a mis-identification of the host galaxy, but 
we show that the large offset is consistent with the merger of 
a NS-NS binary that has been created by dynamical encoun- 
ters within a globular cluster. The offset distribution arising 
from the globular cluster systems of the Milky Way and a 
local massive elliptical galaxy, M87, are consistent with the 
observed offsets of all short GRBs, including GRB 060502B, 
and that the production rate of such binaries may be com- 
parable with that of field DCO binaries. 
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